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ABSTRACT 


One  of  the  major  difficulties  in  attempting  to  apply  known 
queueing  theory  results  to  real  problems  Is  that  almost  always 
these  results  assume  a time  stationary  Poisson  arrival  process, 
whereas  in  practice  the  actual  process  is  almost  invariably 
nonstationary . ^his  paper  we*  considers  single  server  in- 

finite capacity  queueing  models  in  which  the  arrival  process.  (_/^6c>a 
is  a nonstationary  process  with  an  intensity  function  Att)  , i 
tr5^0  , which  is  itself  a random  process  .^..c-We  STippose  that  the 
average  value  of  the  intensity  function  exists  and  is  equal* to 
constant,  cal]#^!*  , with  probability  1. 

4 — We  make confect  ore  the  effect  that  the  **closer  £A(t),t^£)(Q 


some 


is  to  the  stationary  Poisson  process  with  rate  X**then  the  smaller 
is  the  average  customer  delay,  and  then  wo-  verify  the  conjecture  i5 
in  a special  case. 


i/or 


, dt 
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AVERAGE  DELAY  IN  QUEUES  WITH  NONSTATIONARY  POISSON  ARRIVALS 


Sheldon  M.  Ross 


1 . THE  CONJECTURE 


Consider  a single  server.  Infinite  capacity,  queueing  model  In 


which  the  arrival  process  Is  a nonstationary  Poisson  process  with 


Intensity  function  A(t)  , t ^0  ; where  [A(t),  t ^ 0]  Is  Itself  a 


Is  such  that 


exists  and  Is  constant  with  probability  1.  The  above  would  follow,  for 


regeneration  time  T . In  this  case  X would  be  given  by  (see  [2]) 


A(s)ds 


In  addition  we  suppose  that  the  successive  service  times  are  Independently 


chosen  from  some  service  distribution 


The  order  of  service  Is  "first  come  first  served 


Let  us  denote  by  d the  average  amount  of  time  that  a customer  spends 


waiting  In  queue.  We  then  make  the  following  conjecture  concerning  d 


2 


Conjecture: 


d > mn 

- 2(1  - XE[S]) 


It  should  be  noted  that 


would  be  the  value  of  d If  the 


xw  wtiow  2(1  XE[S])  wwuxu  uc  uiic  voxuc  wl  u xx  ufic 

arrival  process  was  a stationary  Poisson  process  with  rate  X . 

In  fact  not  only  do  we  conjecture  the  above,  but  we  also  believe 
that  the  "closer  {A(t),t  ^0}  is  to  the  stationary  process  with  rate 
X " then  the  smaller  d is.  For  example  suppose  that  {A(t),t  ^ 0} 
is  a 2-state  continuous  time  Markov  chain  alternating  between  the  states 
X^  and  X2  and  suppose  that  the  time  It  spends  in  state  X^  during 
each  visit  is  an  exponential  random  variable  with  rate  ca^  , 1 » 1,2  . 
Let  d(c)  denote  the  average  customer  delay  for  this  model.  It  should 
be  noted  that  the  average  arrival  rate  X , given  by 


^1*^2  ^2“l 

01  + 02 


is  Independent  of  the  constant  c which  regulates  how  fast  the  arrival 
rate  changes  between  its  possible  values.  Intuitively  the  larger  c is 
the  "more  stationary"  the  arrival  process  Is  and  thus  we  conjecture  that 


d(c)  is  a decreasing  function  of  c 


lim  d(c) 


2(1  - XEfS]) 
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We  shall  use  Che  usual  notation  F/G/1  to  describe  the  queueing  system 
in  which  the  arrival  process  is  a renewal  process  with  interarrival  dis- 
tribution F and  there  is  a simple  server  having  service  distribution  G . 
The  following  result  was  proven  by  Stoyan  [ 2 ] . 


Proposition  1: 


Consider  the  2 queueing  systems  F/G/1  and  F^/G^/1.  If  F F^  and 

V 

G G^  and  J*xdF(x)  =J*xdF^(x)  then  ^ » where  (®*) 

limiting  distribution  of  delay  in  the  system  F/G/1  (F^/G^/1). 

The  above  can  be  proven  by  showing  Chat  D < where  D and 

n—  n n n 

V 

represent  the  delays  in  queue  of  the  nth  customers  of  the  respective  models. 
This  is  shown  by  induction  on  n , by  using  properties  2 and  3 (which  implies 


that  if  X < Y and  E[X]  ■=  E{Y]  then 


-X  -y| 


and  Che  identity 


D « max  [0,D  + U ] , 
n+1  n n 


where 


U = nth  service  time  - nth  interarrival  time, 
n 

Remark; 

It  is  perhaps  surprising  that  the  analogue  of  Proposition  1 does  not 
hold  when  there  is  more  than  one  server.  To  best  understand  the  forthcoming 
counterexample  first  consider  2 two-server  queueing  systems,  the  servers 
of  the  first  system  having  constant  service  time  of  1 unit,  and  the  servers 
of  the  second  system  having  service  times  that  are  either  0 or  2 - each 
possibility  having  probability  1/2.  Suppose  that  in  both  systems  customers 


I 
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arrive  In  batches  of  size  3,  with  the  times  between  each  batch  being  large. 
Then  it  is  easy  to  see  that  the  average  delay  in  the  first  (constant) 
system  is  1/3  whereas  it  is  1/6  in  the  second  (more  variable)  system. 

The  above  example  yields  a counterexample  even  when  we  require  the 
arrival  process  to  be  a renewal  process.  To  see  this  suppose  that  the 
arrival  process  is  a renewal  process  with  interarrival  times  where 


1° 


with  probability  e 


large  with  probability  1 - e . 


Hence  in  this  renewal  arrival  process  the  customers  will  arrive  in  batches 
of  random  size.  As  batches  of  size  greater  than  3 occur  much  more  in- 

3 

frequently  than  those  of  size  3 [the  probabilities  being  e versus 
2 

c (1  - e)],  it  is  clear  from  the  preceding  argument  that  average  customer 
delay  is  greater  in  the  constant  service  model. 
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3.  PROOF  OF  THE  CONJECTURE  WHEN  X2  = 0 


Consider  the  Infinite  capacity  single  server  model  of  Section  1 where 
A(t)  alternates  between  and  X2  , spending  an  exponential  amount  of 

time  with  rate  ca^  during  each  visit  to  X^  , i « 1,2  . Let  D(c)  denote 
the  limiting  distribution  of  customer  delay  in  this  system,  and  let  d(c) 
denote  the  mean  of  this  distribution. 


Proposition  2: 

If  X2  = 0 then  D(c^)  DCc^)  whenever  i ^2  ’ thus  d(c) 

V 

is  decreasing  in  c . 


Proof : 


It  is  easy  to  see  that  when  X2  = 0 the  arrival  process  is  a renewal 
process.  Let  denote  an  arbitrary  interarrival  time  of  this  renewal 

process.  We'll  prove  the  proposition  by  showing  that  X ^X  when 

‘^l  V ^2 

C2  £ Cj  . To  show  this  let  M denote  an  exponential  random  variable  with 
rate  X,  , and  let  N denote  a random  variable  which,  conditional  with 

M has  a Poisson  distribution  with  mean  c^a^^CM)  . Also  let  , i ^ » 
denote  Independent  (of  themselves  and  of  all  the  other  random  variables 
defined  above)  exponential  random  variables  having  rate  02  . Then  using 
the  notation  X = Y to  mean  that  X and  Y have  the  same  probability 
distribution,  we  have  the  following 


(3.1) 


1 ^1 
X - M + — +•••  + 

^1  *^1  *^1 


The  above  follows  by  interpreting  E as  the  total  amount  of  time  the 


process  spends  in  its  phase  during  an  interarrival  period; 

can  be  interpreted  as  the  number  of  changes  from  X to  X that 


occur  during  this  interarrival  period;  as  V /c  as  the  time  spent  in 


the  X-  phase  during  the  ith  visit 


we  have  that 


where  the  1^  are  independent  (of  each  other  and  of  all  previously  defined 
random  variables)  and  such  that 


1 with  probability  c /c 


0 with  probability  1 


Equation  (3.2)  is  justified  by  the  same  argument  as  used  In  (3.1)  along 
with  the  fact  that  if  each  of  a Poisson  number  of  events  is  independently 
counted  with  probability  p then  the  number  of  counted  events  is  also 


Now  it  is  straightforward  to  check  that 


and  thus  (using  property  2)  we  see  that 


for  all  increasing  convex  f . It  now  follows  from  (3.1)  and  (3.2)  by 


taking  expectation  of  both  sides  of  the  above  inequality,  that 


and  the  result  follows  from  Proposition  l.B 


Before  obtaining  the  limiting  value  of  d(c)  let  us  first  compute 


From  (3.1)  we  have 


the  mean  and  variance  of  the  interarrival  time  X 


Var[X 


Var[X  ] 


We  are  now  ready  to  show  that  as  c ->  ® d(c)  converges  to  its  corre 


spondlng  value  In  the  stationary  case 


Proposition  3 


where  S represents  an  arbitrary  service  time  random  variable 


Proof 


Marshall  [ 3 1 showed  that  for  GI/G/1  queues 


Var[S]  + Var[X  ] + (E[S]  - E[X  ]) 


where  I is  the  length  of  an  idle  period.  Now  an  idle  period  is  a mixture 


of  a random  variable  having  the  same  distribution  as  X and  a random 


and  an  Independent 


variable  having  the  same  distribution  as  the  sum  of  X 


exponential  random  variable  with  rate  co,  (the  probabilities  for  this 


mixture  being  the  respective  probabilities  that  an  idle  period  starts 


11 


In  a or  a X2  phase).  Hence,  It  follows  that 


(3.6) 


I 


< 

St 


X + U 
c 


where  ^ means  stochastically  smaller  than  and  U Is  an  exponential 

St 

% 

random  variable  with  rate  c«2  Independent  of  the  Interarrlval  random 
variable  . Now,  from  the  representation  (3.6)  and  (3.3)  and  (3.4) 
we  see  that 


11m  E[I] 


1 

X 


11m  E[I^]  - 11m  ElX^l 

C->«o  c.-*co  L 


Letting  c approach  Infinity  In  (3.5)  and  substituting  the  above  yields 
the  result. B 

Remark; 

In  this  section  we  have  supposed  that  the  {A(t),t  ^ 0}  process 
alternates  between  the  values  X^  and  X2  ° 0 , spending  an  exponential 
amount  of  time  with  rate  ctt2  In  the  X2  = 0 phase.  However  all  proofs 
would  go  through  In  an  Identical  fashion  If  we  supposed  that  the  time  In 
the  0-phase  had  a distribution  where  G^(x)  = G(cx)  for  some  arbitrary 

distribution  G having  finite  variance.  On  the  other  hand,  however,  It 
Is  unfortunate  that,  to  maintain  renewal  arrivals,  the  assumption  of 
exponential  times  In  the  X^^-phase  Is  essential  for  our  approach. 
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A.  HEURISTICS  AND  FUTURE  WORK 

One  of  the  main  reasons  for  our  belief  In  the  conjecture  of  Section  1 
Is  based  on  the  following  Intuitive  argument:  Suppose  that  [N(t),t  ^ 0] 

Is  a nonhomogeneous  Poisson  Process  with  Intensity  function  X(t)  , 
and  suppose  that,  for  some  positive  constant  X , 


11m 

S-«o 


S 


/ 


X(t)dt 

s 


X . 


In  considering  processes  satisfying  the  above  as  possible  arrival 
processes  for  a waiting  line  model  we  should  first  note  that  for  large  T 
they  all  have  roughly  the  same  distribution  for  the  number  of  arrivals 
by  time  T - namely  Poisson  with  mean  roughly  XT  . The  difference 
arises  in  the  conditional  distribution  of  the  arrival  times.  Specif ically, 
given  that  n arrivals  have  occurred  by  time  T , these  n arrival 
times  have  the  same  joint  distribution  as  do  n Independent  random 
variables  chosen  from  the  common  distribution  F where 


Now  if  we  had  to  choose  a Poisson  (with  mean  XT)  number  of  arrival 
times  according  to  some  distribution  concentrated  over  [0,T]  so  as 
to  minimize  average  customer  waiting  time  then,  intuitively,  it  seems 
reasonable  that  we  would  want  that  distribution  which  does  not  "favor" 


mi. 


ik 
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